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Abstract 

In this article, we take the scalar diquarks as point particles and describe 
them as basic quantum fields, then introduce the SU(3) color gauge inter- 
action and new vacuum condensates to study the nonet scalar mesons as 
■ tetraquark states with the QCD sum rules. Comparing with the conventional 

quark currents, the diquark currents have the outstanding advantage to satisfy 
the two criteria of the QCD sum rules more easily. 
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y£ : 1 Introduction 

>, 

The light flavor scalar mesons present a remarkable exception for the naive quark 
model, and the structures of those mesons have not been unambiguously deter- 
mined yet. The numerous candidates with J PC = ++ below 2 GeV cannot be 
accommodated in one qq nonet, some are supposed to be glueballs, molecular states 
and tetraquark states (or their special superpositions) [H El [3]. The a (980) and 
/o(980) are good candidates for the KK molecular states [I], however, their cousins 
er(600) and k(800) lie considerably higher than the corresponding thresholds, it is 
difficult to identify them as the tttt and ttK molecular states, respectively. There 
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may be different dynamics which dominate the ++ mesons below and above 1 GeV 
respectively, and result in two scalar nonets below 1.7 GeV [HI21E]- The strong at- 
tractions between the diquark states (qq)?; and (qq)s in relative S'-wave may result in 
a nonet tetraquark states manifest below 1 GeV, while the conventional 3 Po QQ nonet 
have masses about (1.2 — 1.6) GeV, and the well established 3 Pi and 3 P2 qq nonets 
with J PC = 1 ++ and 2 ++ respectively lie in the same region. Furthermore, there 
are enough candidates for the 3 P QQ nonet mesons, a (1450), / (1370), iC(1430), 
/o(1500) and/o(1710) 0. 

In the tetraquark scenario, the structures of the nonet scalar mesons in the ideal 
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mixing limit can be symbolically written as [6] 
^(600) = udud, / (980) 

a Q (980) = dsns, a°(980) 
K +(800) = mc/Js, k°(800) 

The four light isospin-| Kn resonances near 800 MeV, known as the k(800) mesons, 
have not been firmly established yet, there are still controversy about their existence 
due to the large width and nearby Kit threshold [5j. 

In general, we may expect constructing the tetraquark currents and studying 
the nonet scalar mesons below 1 GeV as the tetraquark states with the QCD sum 
rules [HIE]. For the conventional mesons and baryons, the "single-pole + continuum 
states" model works well in representing the phenomenological spectral densities, the 
continuum states are usually approximated by the contributions from the asymptotic 
quarks and gluons, the Borel windows are rather large and reliable QCD sum rules 
can be obtained. However, for the light flavor multiquark states, we cannot obtain 
a Borel window to satisfy the two criteria (pole dominance and convergence of the 
operator product expansion) of the QCD sum rules [9]. In Ref.[10j, T. V. Brito et 
al take the quarks as the basic quantum fields, and study the scalar mesons <r(600), 
k(800), /o(980) and a (980) as the diquak-antidiquark states with the QCD sum 
rules, and cannot obtain Borel windows to satisfy the two criteria, and resort to a 
compromise between the two criteria. For the heavy tetraquark states and molecular 
states, the two criteria can be satisfied, but the Borel windows are rather small [llj. 

We can take the colored diquarks as point particles and describe them as the basic 
scalar, pseudoscalar, vector, axial-vector and tensor fields respectively to overcome 
the embarrassment [12]. In this article, we construct the color singlet tetraquark 
currents with the scalar diquark fields, parameterize the nonperturbative effects with 
new vacuum condensates besides the gluon condensate, and perform the standard 
procedure of the QCD sum rules to study the nonet scalar mesons below 1 GeV. 
The QCD sum rules are "new" because the interpolating currents are constructed 
from the basic diquark fields instead of the quark and gluon fields. 

Whether or not the colored diquarks can be taken as basic constituents is of great 
importance, because it provides a new spectroscopy for the mesons and baryons 

mm- 

The article is arranged as follows: we derive the new QCD sum rules for the nonet 
scalar mesons in Sect. 2; in Sect. 3, we present the numerical results and discussions; 
and Sect. 4 is reserved for our conclusions. 



usus + dsds 
usus — dsds 



Qq (980) = usds, 



V2 

udus, k°(800) = usud, k _ (800) = dsud . 

(1) 
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2 The nonet scalar mesons with QCD Sum Rules 



In the following, we write down the interpolating currents for the nonet scalar mesons 
below 1 GeV, 



Jao 



D a (x)D a (x) + U a (x)U a (x) 
D a (x)D a (x) - U a (x)U a {x) 



V2 

J K (x) = S a (x)U a (x), 

J a (x) = S a (x)S a (x), (2) 



where 



U a {x) oc U a {x) = e ahc dl{x)C^s c {x) , 
D a (x) oc D a {x) = e abc ul(x)C l5 s c (x) , 
S a (x) oc S a (x) = t abc ul{x)C^d c {x) , (3) 

the a, b, c are color indices, the C is the charge conjugation matrix, the U a (x), 
D a (x) and S a (x) are basic scalar diquark fields, while the U a (x), D a (x) and S a (x) 
are the corresponding scalar two-quark currents. In this article, we take the isospin 
limit for the u and d quarks, and denote the fields U a (x) and D a (x) as Q a (x). 

For the general color antitriplet bilinear quark-quark fields q(x)q(yjfl and color 
singlet bilinear quark- ant iquark fields q{x)q(y), where the flavor, color and spin 
indexes are not shown explicitly for simplicity, we can project them into a local 
and a nonlocal part, after bosonization, the two parts are translated into a basic 
quantum field and a bound state amplitude, respectively, 

q{x)q{y) -> D (^y^J T D (x - y) , 

q(x)q(y) -> M T M (x - y) , (4) 

where the D(rr) and M(x) denote the diquark and meson fields respectively, the 
ru(x) and Tm(^) denote the corresponding Bethe-Salpeter amplitudes respectively 
[T4"l [15]. In Ref. [16], we study the structures of the pseudoscalar mesons tt, K and 
the scalar diquarks U a , D a , S a in the framework of the coupled rainbow Schwinger- 
Dyson equation and ladder Bethe-Salpeter equation using a confining effective po- 
tential, and observe that the dominant Dirac spinor structure of the Bethe-Salpeter 
amplitudes of the scalar diquarks is C75. If we take the local limit for the nonlocal 
Bethe-Salpeter amplitudes, the dimension-1 scalar diquark fields U a , D a and S a are 



2 If we take the local limit in the scalar channels, the two-quark currents U a (x), D a (x) and 
S a (x) are recovered. 
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proportional to the dimension-3 scalar two-quark currents U a , D a and S a , respec- 
tively. A dimension-1 quantity A can be introduced to represent the hadronization 

The attractive interaction of one-gluon exchange favors formation of the diquarks 
in color antitriplet 3 C , flavor antitriplet 3/ and spin singlet l s [21 [17]. Lattice QCD 
studies of the light flavors indicate that the strong attraction in the scalar diquark 
channels favors the formation of good diquarks, the weaker attraction (the quark- 
quark correlation is rather weak) in the axial-vector diquark channels maybe form 
bad diquarks, the energy gap between the axial-vector and scalar diquarks is about 
| of the A-nucleon mass splitting, i.e. ~ 0.2 GeV [18], which is also expected from 
the hypersplitting color-spin interaction m l m Ti ■ TjBi ■ <7j [SI HZ]- On the other hand, 
the studies based on the random instanton liquid model indicate that the instanton 
induced quark-quark interactions are weakly repulsive in the vector and axial- vector 
channels, strongly repulsive in the pseudoscalar channel, and strongly attractive in 
the scalar and tensor channels [19]. So it is sensible to use the scalar diquark fields 
to construct the tetraquark currents. 

The two-point correlation functions Ilj(p) can be written as 

Ik{p) = i [d 4 x e^(0|T[J,(x)J, t (0)]|0) , (5) 



where the current Ji(x) denotes Jf (x), J ao (x), J K (x) and J a (x). 

We can insert a complete set of intermediate hadronic states with the same 
quantum numbers as the current operators Ji(x) into the correlation functions Hi(p) 
to obtain the hadronic representation [7J [8] . Isolating the ground state contributions 
from the pole terms of the nonet scalar mesons, we obtain the results, 

rr M Aj ^ (0|^0)l%))(%M(0)|0) 

l[P) Mf-p 2 ^^ M 2 -p 2 



.702 S~p 2 J s S-p 2 

where the Mj are the ground state masses, the Aj are corresponding pole residues de- 
fined by (0\Ji(0)\S(p)) = X h the thresholds 6^ = (2m s + 2m g ) 2 , (m s + 3m q ) 2 , (Am q ) 2 
in the channels /o/ a o(980), k(800), a(600) respectively, the s° are the thresholds for 
the higher resonances and continuum states \h), and the p^{s) are the corresponding 
hadronic spectral densities. 

We introduce the following new Lagrangian £, 

2 2 M 2 M 2 u 2 D 

- 1 -m^S + ..., (7) 
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Figure 1: The typical diagrams we calculate in the operator product expansion, 
where the solid line denotes the di quark propagators. 

where = <9 M — ig s G^ and carry out the operator product expansion with the 
basic diquark fields U, D and S instead of the quark fields u, d and s, and we have 
neglected the terms concerning the heavy diquark fields in the Lagrangian. In the 
QCD, the basic quantum fields are the quark and gluon fields, the attractive interac- 
tions in the color antitriplet 3 C , flavor antitriplet 3/ and spin singlet l s quark-quark 
channels favor formation of the scalar diquarks, we can absorb some effects of the 
quark-gluon interactions into the effective diquark masses, which are characterized 
by the correlation length L ~ ^. At the distance / > L, the 3 C diquark state com- 
bines with one quark or one 3 C antidiquark to form a baryon state or a tetraquark 
state, while at the distance / < L, the 3 C diquark states dissociate into asymptotic 
quarks and gluons gradually, the strength of the quark-quark correlations is very 
weak. Just like the quarks, the diquarks have three colors, and can be gauged with 
the same SU (3) color group to embody the residual quark-gluon interactions. 

In calculations, we take into account all diagrams like the typical ones shown 
in Fig.l, introduce new vacuum diquark condensates (QQ) and (SS) besides the 
gluon condensate to parameterize the nonperturbative QCD vacuum, and consider 
the vacuum condensates up to dimension four. In the QCD sum rules, the high 
dimensional vacuum condensates are usually suppressed by large denominators or 
additional powers of the inverse Borel parameters the net contributions are 
very small. For example, in the present case the contributions of the dimension-6 
vacuum condensates can be estimated as {QQ){ asGG )^ ~ aoooo ^ GcV ; which is a 
tiny quantity. If additional suppressions originate from the denominators are taken 
into account, the contributions are even smaller, and can be safely neglected. 
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Once the analytical results are obtained, then we can take the dualities below 
the thresholds Sq and perform the Borel transform with respect to the variable 
P 2 = —p 2 , finally we obtain the following sum rules, 

\ 2 e~^ = / dse~^pi(s), (8) 

where the Pi(s) denote the QCD spectral densities Pf / ao (s), Po-( s ) an d Pk{s), 

I \ ^ n//S^\ c-/ 2 \ ™Q ,Ot s GG If-/ ~2 \ 

Pfo/aois) = + 2<QQ>5(* - m Q ) - Qttfii—^) J o dx ^{s ~ m Q ) , 

PM = ^ + 2(SS)5(s-m*)-^(^) ^ dx- 3 S(s - m 2 ) , 

lb'K z 96i 4 7T Jo x 



p K (s) = -^ + (QQ)6(s~m 2 s ) + (SS}6(s-m 2 Q 



lGir 



1 ,a*GG s f 1 , \m 2 Q m 2 s 



192T 4 7T 



) / dx 



o 



+ 



x 3 (1 — x) 3 



S{s - m% /s ) , (9) 



m s/Q = ^f)' ^s/Q = -^ + S> A ? = ( 2m Q) 2 ' ( m Q + m 5) 2 , (2m s ) 2 in the channels 
/o/ao(980), k(800), cr(600) respectively, and the T 2 is the Borel parameter. The 
threshold parameters A 2 are different from the corresponding 0^ in Eq.(6), because 
we absorb some QCD interactions into the effective diquark masses. 

Differentiate Eq.(8) with respect to then eliminate the pole residues Aj, we 
can obtain the sum rules for the masses of the nonet scalar mesons, 

o s 

m 2 = ff a ^FW^> e T (10) 

J^ 2 dspi(s)e 5* 



3 Numerical Results 

We estimate the vacuum diquark condensates (QQ) and (55) with the assumption 
of the vacuum saturation, which works well in the large N c limit, 

(QQ) = { -^ = AHQQ), 

(SS) = ^^ = A 4 (55), (11) 
6 

where the A is a quantity has dimension of mass and can be taken as the confinement 
energy scale A = 0.3 GeV. At the energy scale p = 1 GeV, (ss) = 0.8(qq), (qq) = 
-(0.24 GeV) 3 , = (0.33 GeV) 4 [Zl IS], and (QQ) = 0.8(55) = 0.0031 GeV 2 . 

The quark condensates play a special role being responsible for the spontaneous 
breaking of the chiral symmetry, and relate with the masses and decay constants of 
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the light pseudo-scalar mesons through the Gell-Mann-Oakes-Renner relation [20J. 
The values of other vacuum condensates, such as the mixed condensates, the four 
quark condensates and the gluon condensates, cannot be obtained from the first 
principles, we usually calculate them with the lattice QCD, the instanton models, or 
determine them empirically by fitting the QCD sum rules to the experimental data. 
In this article, we introduce the diquark condensates (QQ) and (SS) to parameterize 
the nonperturbative QCD vacuum, and assume that they relate with the four quark 
condensates (therefore they have implicit relations with the spontaneous breaking of 
the chiral symmetry), and take the dimension-one parameter A to be the confinement 
energy scale, as the scalar mesons /o(980), ao(980), k(800) and a(600) are bound 
states which consist of the confined quarks and gluons, and some parameters are 
needed to embody the confinement. On the other hand, we can understand the 
parameter A = 0.3 GeV as a fitted value, which happens to be the confinement 
energy scale, the crude estimation works well. 

We take the updated values of the diquark^ masses^ from the QCD sum rules for 
consistency, where the interpolating currents U a (x), D a (x) and S a (x) are used [2Tj . 
rriQ = 0.46 GeV and m s = 0.40 GeV; the scalar diquarks were originally studied with 
the QCD sum rules about twenty years ago [22]. There have been several theoretical 
approaches to estimate the diquark masses, for example, the simple constituent 
diquark mass plus hyperfine spin-spin interaction model [23J. 

The /o(980) and ao(980) are well established, and the existence of the cx(600) me- 
son is confirmed, although there are controversy about its mass and width, the values 
listed in the Review of Particle Physics are (400 - 1200) MeV and (600 - 1000) MeV 
respectively [5J. As far as the k(800) are concerned, there are still controversy 
about their existence, we take them as the 5*-wave isospin-| Kn resonance with 
the Breit-Wigner mass about 850 MeV. The E791 collaboration observed a low- 
mass scalar Kit resonance with the Breit-Wigner mass M = (797± 19 ± 43) MeV 
and width T = (410 ± 43 ± 87) MeV respectively in the decay D + -> K-n + n + 
[24], and the BES collaboration observed a clear low mass enhancement in the 
invariant Kir mass distribution in the decay J ftp — > K*(892)K + ir~ with the Breit- 
Wigner mass M = (878 ± 23^) MeV and width T = (499 ± 52tg) MeV, respec- 
tively [25]. Recently, the BES collaboration reported the charged k(800) in the 
decay J/ip -> K*(892)^K s tt ± with the Breit-Wigner mass M = (826 ± 49l£t) MeV 
and width T = (449 ± 156^gf 4 ) MeV, respectively [26]. It is sensible to estimate 
Mf / ao — M K = m s — m q = 0.14 GeV. On the other hand, the QCD sum rules for the 
tetraquark states indicate that M K = (0.80-0.88) GeV and M a = (0.72-0.80) GeV 

Assuming the energy gap between the ground and first radial excited tetraquark 
states is about 0.5 GeV, we can tentatively determine the threshold parameters 
s h/ao = (i.o + 0.5) 2 GeV 2 , s° K = (0.85 + 0.5) 2 GeV 2 , and s° = (0.75 + 0.5) 2 GeV 2 . 

The convergence behavior of the operator product expansion is very good, the 
contributions from the different terms have the hierarchy: perturbative-term > 
(QQ) S> ( a )■ In calculation, we take uniform minimum value for the Borel 
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parameters T^ in = jj? = 1.0 GeV 2 . The perturbative continuum y^2@( s — s o) * s 

s 

suppressed by the factor e t 7 , the contributions from the pole terms are very large, 
see Fig.2. In this article, we take uniform maximum value for the Borel parame- 
ters, T^ ax = 1.9 GeV 2 , the contributions from the pole terms are about (61 — 84)%, 
(54 - 79)% and (51 - 75)% in the channels / (980)/a (980), k(800) and (x(600), 
respectively. The two criteria (pole dominance and convergence of the operator 
product expansion) of the QCD sum rules are well satisfied. 

If the conventional quark currents are chosen, the multiquark states i.e. tetraquark 
states, pentaquark states, hexaquark states, etc, have the spectral densities p ~ s n 
with the largest n > 4, the integral f Q s n e~T*ds converges slowly [9]. If one don't 
want to release the criterion of pole dominance, we have to either postpone the 
threshold parameter sq to very large value or choose very small value for the Borel 
parameter T^ ax . With large value of the threshold parameter s , for example, 
s ^> M 2 r , here gr stands for the ground state, the contributions from the high 
resonance states and continuum states are included in, we cannot use single-pole (or 
ground state) approximation for the spectral densities; on the other hand, with very 
small value of the Borel parameter T^ ax , the Borel window T^ ax — T^ in shrinks to 
zero or very small values. 

The numerical values of the masses and pole residues are presented in Table 1 
and Figs. 3-4. From Table 1, we can see that the present predictions are compatible 
with (or not in conflict with) the experimental data [5l [2H [251 12H] and theoretical 
estimations [27J. 

The scalar tetraquark currents Jf / ao ( x )-> Jk( x ) and J c {x) maybe have non- 
vanishing couplings with the scattering states tttt, KK, Kn, Kr], ttt], rjr], etc, for 
example, 

(0|J /o (0)|7T7r(p)) = X foWW , 
(0\J fo (0)\KK(p)) = X foKK . (12) 

If the couplings denoted by the A f on7T and A f KK are strong enough, the contamina- 
tions from the continuum states are expected to be large. In the following, we study 
the contributions of the intermediate pseudoscalar meson loops to the correlation 
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function Hf (p) in details as an example 



1 2 _ f ° M 2 9foKK^KKip)g fo KK * 
^ 70 r 70 



P _M /o 



p2 


- M? 






^ 2 


-Ml 






p2 


-Ml 

/() 






^2 





gfoTrTr^Tr(j>)T ww ^ K K^KK(P)9foKK-^—J72- 

■P fo 



9foKK^KK(p)T K K^mT^nn(p)gfo 



^fo 



7777 P 2 ~ M% 

9f KK^KK(p)T K K^KK^KK(p)9foKK 2 _ , , 2 '(^) 

P iW /o 



where 



:(P) 



(27r) 4 [q 2 — mj] [(p — q) 2 — m 
d 4 q 1 



21 



^ KK{P) J W[q 2 -m 2 K ][{p-q) 2 -m 2 K Y {U) 

the g/ 07rvr , gf KK are the strong coupling constants between the / (980) and the 
pseudoscalar meson pairs tctt, KK respectively, and the T nn ^ nn , T^-^kRi T^kk-^-k-k-, 
Tkk^kk are the scattering amplitudes among the pseudoscalar meson pairs ixix 
and KK. The couplings A/ 07rvr and X/ kk are complicated functions of the M/ , p 2 , 
A/,,, Tnn-ymr, T nT ^ KR , T KR ^ nw , T kR _+ kR , g fomr and g fo K K , the explicit expressions 
are difficult to obtain. We should bear in mind that the intermediate meson loops 
contribute a self-energy to the scalar meson / (980), and therefore the scalar me- 
son / (980) develops a Breit-Wigner width. In fact, the scalar mesons / /°o(980), 
k(800) and cr(600) below 1 GeV can be generated dynamically from the unitaried 
scattering amplitudes of the pseudoscalar mesons [28]. We can take into account 
those meson loops effectively by taking the following replacement for the hadronic 
spectral density, 

Sis -Mi) -> ^# , (15) 

here we neglect the complicated renormalization procedure, take the physical values, 
and ignore the energy scale dependence of the mass Mf and pole residue A /„ for sim- 
plicity; the approximation works well. The QCD sum rules for other scalar mesons 
are treated with the same routine. In Ref.[29], Dai et al perform the renormaliza- 
tion procedure in details to take into account the contributions from the continuum 
states. 
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tetraquark states 


mass (GeV) 


pole residue (GeV) 


/ (980)/a (980) 


c\ r\ A -i r\ A 

0.94 — 1.04 


0.16 - 0.17 


«(800) 


0.84-0.93 


0.15-0.16 


a(600) 


0.77-0.84 


w 0.15 


/o(980)/a (980) [(40 - 100) MeV] 


0.94 - 1.04* 




k(800) [550 MeV] 


0.79 - 0.88* 




a(600) [400 MeV] 


0.77 - 0.84* 




a(600) [800 MeV] 


0.72 - 0.80* 




<x(600) [1200 MeV] 


0.69 - 0.77* 





Table 1: The masses and the pole residues of the nonet ++ tetraquark states. The 
star denotes the modified masses from the sum rules where the finite widths shown 
in the bracket are taken into account in the hadronic spectral densities. 

The widths listed in the Review of Particle Physics are Tf — (40 — 100) MeV, 
r ao = (50-100) MeV, T K = (550±34) MeV and T a = (400-1200) MeV, respectively 
[5J. Taking into account the finite widths, we can obtain the modified masses from 
the corresponding QCD sum rules, which are shown in Table 1 and Fig.5. In cal- 
culation, we observe that the narrow width r/ / ao modifies the mass Mf / ao slightly 
and the effect can be neglected safely, while the broad widths T K and r CT reduce the 
masses M K and M a about 55 MeV and (0 — 75) MeV, respectively. Comparing with 
the experimental data [HI EH EH1 [26] , the modified masses are better. 

The scalar tetraquark currents Jf / ao (x), Jk(%) an d Ja{x) maybe also have non- 
vanishing couplings with the higher resonances ao(1450), /o(1370), if* (1430), /o(1500) 
and /o(1710), which are supposed to be the 3 Pq qq states or glueballs, the couplings 
should be very small. Furthermore, the threshold parameters are so = 2.25 GeV 2 , 
1.82 GeV 2 and 1.56 GeV 2 in the channels / /a (980), k(800) and a(600), respec- 
tively, the contaminations should be very small. 

4 Conclusions 

In this article, we take the scalar diquarks as the point particles and describe them 
as the basic quantum fields, then introduce the SU(3) color gauge interactions, 
and construct the tetraquark currents which consist of the scalar fields to study 
the nonet scalar mesons as tetraquark states with the new QCD sum rules. The 
numerical values are compatible with (or not in conflict with) the experimental data 
and theoretical estimations. Comparing with the conventional quark currents, the 
diquark currents have the outstanding advantage to satisfy the two criteria of the 
QCD sum rules more easily, the new sum rules can be extended to study other 
multiquark states. 
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Figure 2: The contributions from the pole terms. The A, B and C correspond to 
the channels / (980)/a (980), k(800) and a(600) respectively. 
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Figure 3: The masses of the nonet scalar tetraquark states. The A, B and C 
correspond to the channels /o(980)/a (980), k(800) and <r(600) respectively. 
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Figure 4: The pole residues of the nonet scalar tetraquark states. The A, B and C 
correspond to the channels /o(980)/ao(980), «(800) and <t(600) respectively. 



1.6 
1.5 
1.4 
1.3 
> 1.2 

CD 1.1 



A [40-100 MeV]; 

B [550 MeV]; 

C [400 MeV]; 

C [800 MeV]; 

C [1200 MeV]. 




Figure 5: The modified masses of the nonet scalar tetraquark states. The A, B and 
C correspond to the channels /o(980)/a (980), k(800) and cr(600) respectively. The 
values in the bracket denote the finite widths in the hadronic spectral densities. 
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